Self consistent determination of plasmonic resonances in ternary nanocomposites 
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We have developed a self consistent technique to predict the behavior of plasmon resonances in 
multi-component systems as a function of wavelength. This approach, based on the tight lower 
bounds of the Bergman-Milton formulation, is able to predict experimental optical data, including 
the positions, shifts and shapes of plasmonic peaks in ternary nanocomposites without using any 
fitting parameters. Our approach is based on viewing the mixing of 3 components as the mixing of 2 
binary mixtures, each in the same host. We obtained excellent predictions of the experimental optical 
behavior for mixtures of Ag:Cu:Si02 and alloys of Au-Cu:SiC>2 and Ag-Au:H20, suggesting that the 
essential physics of plasmonic behavior is captured by this approach. 



Accurate control and prediction of plasmonic behavior in metal nanoparticles in various configurations 
are expected to realize future ultra-high density integrated photonic devices [ly. For example, plasmonic ef- 
fects in metal nanoparticle waveguide arrays and metal-dielectric interfaces have been used to transport 
electromagnetic energy below diffraction limited dimensions Glasses containing metal 

nanoparticles exhibit intensity dependent refractive index values several orders of magnitude greater than 
that of silica glass due to dielectric and quantum confinement effects [7]. These enhanced effects, occurring 
due to surface plasmon effects at the interface of the metal and the dielectric, could allow various use of 



inear and nonlinear optical responses of mixtures 



provided appropriate materials are selected 



11]. Experimental observations have also demonstrated that composites containing multiple metals show 



multiple plasmonic peaks as well as dependence of such peaks on metal fraction 1 12J, 13|]. Therefore, the 
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large optical nonlinearities, the waveguiding capability and multiple-wavelength sensitivity offered by plas- 
mon behavior makes it an important area of study with potential applications to optoelectronics, sensing, 
etc. In this Letter, we develop and demonstrate a self consistent technique to accurately predict the effective 
dielectric constant of multi-component systems and apply it to ternary nanocomposites for which reliable 



experimental data are available [1 



3Q 



The accurate prediction of the optical properties of nanocomposites made from multiple optical mate- 
rials is an outstanding problem in current research. This is primarily because the effective permittivity is 
not uniquely determined by the optical properties of the individual components but also requires knowledge 



of the composite microstructure I15L 



16]. Previous theoretical work in this area has focused primarily on 



estimating the effective permittivity of two-component mixtures where one of the components is a dielectric 



and the other is a metal 



ponent composites ll8| 



19L 




. There are numerous models that predict the effective permittivity of two com- 
20|]. Moreover, there are already established rigorous bound for the values of the 



permittivity in such composites 



21, 



2211 . However, these models have not been extended to or applied 



towards interpreting the behavior of multi-metal composites and typically multi-parameter fitting routines 



have been used to reproduce experimental data 



23Q. But such fitting approaches are of limited use in 



predicting the fundamental behavior of unexplored nanocomposite systems and hence, are not generally 
applicable tools for new product design. Recently it was shown that the problem of calculating the bounds 
for a multicomponent system can be reduced to the determination of the convex-hull generated in the com- 



plex plane when the volume fraction of each component act as a barycentric coordinate for the hull [24]. 
One of the limitations of this approach is that it can be used only for obtaining the loosest bounds which 
are the Wiener bounds [25]. These loose bounds cannot predict the location of the plasmon resonances in 
frequency space for a simple two component system. This is because the resonance enhancement of the 
absorption takes place at / m (e^ — ei) + ei = and is in contradiction to the well known condition for the 
resonance at 2e/i + e\ = 0, where f rn is the volume fraction of the metal, eh is the dielectric permittivity of 
the host and e\ is the real part of the dielectric permittivity of the metal or inclusion. 

Here we report a self consistent technique that accurately predicts the plasmonic behavior, including 
peak positions, shifts and shapes, in multi-metal mixtures and/or alloy nanocomposites without the need 
for any fitting parameters. Our approach begins with the fact that for a two component nanocomposite the 
tightest lower bound on the dielectric permittivity can be expressed through the Bergman-Milton formula 



[15 



16], (which is based on restricting the values of the effective dielectric permittivity in the complex 



plane using suitably defined conformal mappings) as: 



(7) 
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where f a and fh denote the volume fractions of the constituent materials with relative permittivity e a and 
eh respectively, and f a + fh = 1- The parameter 7 takes the values 2(1 — / a )/3 < 7 < 1 — 2/ a /3. In this 
work we use 7 = 1 — 2/ a /3 because it corresponds to the exact Maxwell Garnett mixing rule [26]. The 
upper bound corresponds to the complementary mixture and is not of interest to us. 

The central idea of this paper is to use Eq. [TJ in a self consistent fashion within the effective medium 
approximation to calculate the effective permittivity of a three component or ternary nanocomposite. We 
begin with the hypothesis that a three component system formed by e a , e^ and the host can be viewed as 
a mixture of two components, each having an effective permittivity e"^ and calculated using Eq. [TJ 
In this mixing process the average electric field within the composite is held fixed at every stage of mixing. 
The constraint in the volume fraction is introduced while calculating the individual effective permittivities 
while the final effective permittivity is calculated using equal volumes of each mixture. Moreover, notice 
that eh, which represents the dielectric permittivity of the host matrix, is common to both mixtures. This 
simple binary mixing rule, depicted graphically in Fig. [TJ can be stated as follows: the effective permittivity 
of an N -component mixture can be determined by mixing N-l binary mixtures, each comprising of a host 
and a distinct metal, with the host being common to the N-l pairs. In this Letter we focus on applying this 
mixing rule to ternary systems for which experimental data exists. 

In order to test the validity of our hypothesis we have calculated the optical absorption coefficient of 
ternary composites containing multi-metal and/or alloy nanocomposites using experimental data for the 
individual permittivities of the metals obtained from ref. |27|. The binary mixing rule was applied to the 
formation of nanocomposites made from mixtures of two metals, as shown in Fig. [Ha), and for alloys, as 
shown in Fig. |TJb). We also introduced a modification to the imaginary component of the permittivity of 
each metal in order to account for the enhanced rate of electron scattering due to particle size-dependent 



effects by modifying the Drude model 1280. This was accomplished by expressing the imaginary component 
in the high frequency limit as follows: 

3 / \ 

U P _ Jrulk ( 2d + V FT bulk \ 



w 2 r e// \ 2d J 

where uj p is the bulk plasmon resonance frequency of the metal, and r e ff is an effective relaxation time 
given by: 

1 1 + £ (3) 



T eff Uulk 2d 

where r&^fc is the bulk relaxation time of the electron, vp is the speed of the electrons close to the Fermi 
surface, and d is the nanoparticle diameter. The second term in Eq. [3] takes into account the collision rate 
of the electron with the metal walls for spherically shaped particles. 
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In order to compare our predictions to previous experimental measurements, we identified multi-metal 
composites for which experimental observations of the plasmonic behavior are available including the three 
parameters required for the calculations, namely: (i) the average metal particle size; (ii) the length of the 
sample which determines the absorption, and (iii) the composition of the metal particles. Based on this, we 
compared theoretical predictions with the experiments of Magruder and co-workers on Ag and Cu mixtures 
in Si0 2 H by Battaglin and co-workers on Au-Cu alloys in Si02 [13] and by Moskovits and co-workers 



on colloidal solutions of Au-Ag alloy nanoparticles in H2O 11411 . In our calculations, the optical absorption 
in a silica glass or water matrix was estimated in terms of the optical density (OD) given by OD = al/2.3, 
where a is the absorption coefficient of the final mixture given by a = ^ J ^Im(e e ff(j)), where n Q is the 
dielectric constant of the medium (corrected for addition of the nanoparticles), A is the incident wavelength 
and and / is the length of the sample. 



Magruder et al. II12H measured the behavior of the plasmonic peaks as a function of the volume fractions 
of the individual metals in a composite of Ag and Cu in Si02- Here, the various ratios of the metal ions 
were achieved by sequential ion implantation. The resulting composites had particles with a mean size of 
d = 30 nm for Ag and Cu in ratios of 9:3, 6:6, and 3:9 respectively. The central observations were that: (i) 
a dominant peak, attributed to Ag, shifts from ~410 nm to 440 nm as the Ag concentration was increased; 
and (ii) a second weak peak near ~575 nm appeared as the Cu concentration increased. Using the above 
values for the particle diameter and individual metal concentrations along with the experimental values for 
vp and Tbuik for Cu and Ag we calculated the effective permittivity using the mixing rule developed above 
and estimated the plasmonic behavior of this system. Fig. [2]compares the experimental data extracted from 
ref. [12] (symbols) with the results of our predictions (lines). The theory, without the use of any fitting 
parameters, predicts the position of the plasmon peaks as well as the shift occurring due to change in Ag 
concentration. Our calculations predicted a shift from 417 nm to 452 nm with increasing Ag concentration. 
The calculation also accurately captured the appearance of the second peak at 572 nm with increasing Cu 
concentration. Moreover, the widths of the individual peaks, which is extremely sensitive to particle size, 
were fairly well represented by our model. The peak widths were a consequence of the correction introduced 
to the electron scattering time based on the particle size. 

In Fig. [3] we compare model predictions with the data of Battaglin et al. I13I in which the plasmonic 
behavior was studied as a function of alloy concentration and particle size for Au-Cu alloys in Si02- Equal 
ratios of Au and Cu were implanted and then subjected to different annealing conditions. The as-deposited 
sample contained an alloy with AU2CU1 composition (with the remaining Cu in atomic state) and average 
particle size of 3.8 ± 1.2 nm. The sample annealed in H2 had a AuiCui alloy with an average particle 
size of 8.7 ± 2.5 nm while the sample annealed in air contained primarily Au nanoparticles of 33 ± 15 nm 
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diameter, with the Cu preferentially found as an oxide in the near surface of the samples. For the theoretical 
calculations we first created the appropriate alloys using the mixing rule (depicted in Fig. [lib)) and then 
the final composite based on the volume fraction of the alloy and the experimentally assigned particle 
diameters. The results are shown in Fig. [3l with experimental data as symbols and theory as lines. Once 
again the theory accurately predicts the peak positions and shift for this alloy system without the need for 
any fitting parameters. 

In Fig. |4] we compare the theoretically predicted and experimentally observed behavior of colloidal 
solutions containing Au-Ag core-shell nanoparticle alloys of various compositions in a water matrix, as 
measured by Moskovits and co-workers Jl4ll . Their primary conclusions were: (i) a single plasmon peak 
appeared for the core-shell structures and this peak shifted with changing Au-Ag concentration. This behav- 
ior was attributed to the formation of an alloy in the shell whose composition changed with varying Au-Ag 
fraction; and (ii) the trend could be fitted by applying a multi-parameter fitting routine to obtain the best fit 
based on summation of Lorentzian peaks. We applied our mixing rule to first create the appropriate alloy 
(Fig. QIb)) with the particle diameters taken from 11411 . The resulting theoretical behavior is shown as lines 
in Fig. U] The theory again predicted quite well the shift in the position of the plasmon peak with varying 
alloy composition. 

In conclusion, we have developed a self consistent technique to determine/predict the effective permit- 
tivity of ternary composites containing mixtures and/or alloy nanocomposites using the Milton-Bergman 
lower bound expression for two component composites. We have successfully tested our theory against 
previously studied experimental systems comprising mixtures of multi-metal nanoparticles and/or alloys 
of Ag:Au:Si02, Ag-Cu:Si02, and Au-Cu:H20. This approach predicts quite accurately the peak position 
and shift of plasmonic behavior in these mixtures and alloys by utilizing experimentally available optical 



parameters 112711 and without the need for any fitting. This result clearly indicates that our mixing approach 
captures the essential physics of plasmonic behavior in these multi-component systems. This calculation 
can readily be extended for nanocomposites with more than three components. However, at this point, no 
experimental data exists to validate the calculations for such systems. In addition, the shape of the plasmon 
peaks could also be fairly well reproduced by utilizing a simple correction to the relaxation time of the elec- 
trons that accounts for the enhanced scattering at the metal boundary as the particle size is reduced. This 
model allows for further corrections due to quantum confinement effects as the particle size gets reduced 
below the 10 nm size scale |29[]. This approach could be of great value towards predicting optical proper- 
ties in terms of plasmonic behavior as a function of volume fraction, particle size and alloy composition in 
multi-component mixtures and could guide the assembly of nanocomposites with tailored optical properties 
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Figure captions 

1. Schematic of the binary tree approach to implement the self consistent binary mixing rule to obtain 
the effective dielectric constant, (a) Application of the mixing rule to create a mixture of two metals 
a and b in the host h. (b) Application of the rule to create an alloy nanocomposite in host h from 
alloys of metals a and b. 

2. Comparison of theory and experiment for plasmonic behavior in a nanocomposite containing various 
mixtures of Au and Cu in Si02- The theoretical results using the binary mixing rule are presented as 



lines while experimental data (symbols) were extracted from the work of Magruder et al. [12]. The 
peak position, shift and widths are predicted very well using the mixing rule. The experimental data 
corresponds to various ratios of Ag:Cu achieved by sequential ion implantation. 

3. Comparison of theory and experiment of plasmonic behavior in various Au-Cu alloy nanocomposites 
in Si02 obtained by ion implantation and annealing. The experimental data (symbols) was extracted 



from the work of Battaglin et al. H 1 3M while the theoretical calculations are represented as lines. 



Alloy compositions are indicated on the figure. 

4. Comparison of plasmonic behavior in colloidal solutions containing Au-Ag core-shell alloys for two 
compositions. The symbols are experimental data extracted from Moskovits et al. J14I1 while lines 
are theory. The numbers within brackets in the legends represent the mole fraction of Au in the 
colloidal solution. 
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(a) (b) 



Figure 1 : Schematic of the binary tree approach to implement the self consistent binary mixing rule to obtain the 
effective dielectric constant, (a) Application of the mixing rule to create a mixture of two metals a and b in the host h. 
(b) Application of the rule to create an alloy nanocomposite in host h from alloys of metals a and b. 
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Figure 2: Comparison of theory and experiment for plasmonic behavior in a nanocomposite containing various mix- 
tures of Au and Cu in SiCV The theoretical results using the binary mixing rule are presented as lines while experi- 
mental data (symbols) were extracted from the work of Magruder et al. IU2I1 . The peak position, shift and widths are 
predicted very well using the mixing rule. The experimental data corresponds to various ratios of Ag:Cu achieved by 
sequential ion implantation. 
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Figure 3: Comparison of theory and experiment of plasmonic behavior in various Au-Cu alloy nanocomposites in 
Si02 obtained by ion implantation and annealing. The experimental data (symbols) was extracted from the work of 
Battaglin et al. Ill 311 while the theoretical calculations are represented as lines. Alloy compositions are indicated on 
the figure. 
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Figure 4: Comparison of plasmonic behavior in colloidal solutions containing Au-Ag core-shell alloys for two compo- 
sitions. The symbols are experimental data extracted from Moskovits et al. IU4I1 while lines are theory. The numbers 
within brackets in the legends represent the mole fraction of Au in the colloidal solution. 



